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Abstract. In this work, we describe a possible experimental realization of Bose's 
idea to use spin chains for short distance quantum communication [S. Bose, Phys. Rev. 
Lett. 91 207901]. Josephson arrays have been proposed and analyzed as transmission 
channels for systems of superconducting charge qubits. Here, we consider a chain 
of persistent current qubits, that is appropriate for state transfer with high fidelity 
in systems containing flux qubits. We calculate the fidelity of state transfer for this 
system. In general, the Hamiltonian of this system is not of XXZ-type, and we analyze 
the magnitude and the effect of the terms that do not conserve the z-component of the 
total spin. 
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1. Introduction 

Recently, the idea to use quantum spin chains for short-distance quantum 
communication was put forward by Bose lj. He showed that an array of spins (or 
spin-like two level systems) with isotropic Heisenberg interaction is suitable for quantum 
state transfer. The advantage of spin chains as transmission lines is the fact, that they 
do not need to have controllable couplings between the qubits or complicated gating 
schemes to achieve high transfer fidelity. An initial state is prepared at one end of the 
chain at time t = 0, and after a certain time ti is measured at the other. The fidelity 
of quantum communication averaged over all pure input states on the Bloch sphere is 
taken as a measure of the transmission quality. 

Bose showed that for short chains (number of spins ~ 100) the average fidelity is 
quite high, greater than 2/3, which is the highest fidelity of transmission through a 
classical channel j2J. In a homogeneous chain, i.e. if all coupling constants are the same, 
the information about the input state is dispersed between the spins at all times t > 0. 
Therefore the fidelity is always less then unity. 

Some methods were proposed to achieve perfect state transfer with fidelity one. 
A special form of the Hamiltonian with spatially varying coupling constants between 
the qubits allows to avoid dispersion [3J H] . Another method is to form Gaussian wave 
packets (with low dispersion) by encoding the information using multiple spins . Also 
the combination of two spin chains jH] can be used to achieve perfect state transfer. This 
method has the advantage that it can be implemented using almost any two spin chains 
and is stable to fluctuations of the chain parameters |7|. However, the time after which 
perfect state transfer is achieved grows if the individual fidelities of the chains decrease. 
It is therefore advantageous to have single chains with high fidelity to implement this 
improved method. 

Quantum state transfer can be implemented using any type of two-level systems. 
However, it is preferable to use a technology that is adapted to the quantum information 
hardware that is supposed to be coupled by the transmission line. One of the most 
promising architectures of quantum computing devices are superconducting circuits, 
for example charge, flux and charge-flux qubits. In recent years these were intensively 
studied both theoretically and experimentally. 

A possible realization of an effective transmission line for charge qubits was 
described in jH]. There, the fidelity of state transfer through Josephson junction arrays 
and the influence of static disorder and dynamical noise was analyzed. 

2. Arrays of persistent-current qubits 

In this work we consider a line of persistent-current qubits [9]. We will show that it 
is appropriate for state transfer with high fidelity in systems containing flux qubits. A 
persistent current qubit [TU] is a superconducting loop with three Josephson junctions, 
see Fig. We assume that the left and right Josephson junctions have capacitance 
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Figure 1. Persistent-current qubit and Josephson energy as a function of the relative 
phase of left and right Josephson junctions. 




Figure 2. Capacitively coupled qubits. 



C and Josephson energy Ej, the central junction is characterized by a capacitance aC 
and Josephson energy aEj with a < 1. The gate capacitances (not shown in the figure) 
are equal 7C. The Hamiltonian of the qubit 

H = -A a x - Ba z (1) 

is the same as that of a spin-| particle in a magnetic field. The eigenstates |0) = | |) 
and |1) = I |) of a z correspond to clockwise and counterclockwise currents. The 
coefficient A is a tunneling amplitude between these states and B depends on the 
flux through the qubit $ and the modulus of the circulating current I p 

B = /„($) ($ - i$ ) , (2) 

here, $ — h/(2e) is the flux quantum. The circulating current I p depends on the 
magnetic frustration, i.e. the amount of external magnetic flux in the loop in units 
of the flux quantum [TI]. The effective magnetic field B is determined by the qubit 
parameters and the external magnetic flux. 

We assume, that the temperature is low enough, i.e. ksT is smaller than the energy 
of the state |1), so we can neglect thermal fluctuations. 

Persistent-current qubits can be capacitively coupled (with coupling capacitance 
(3C, see Fig. |2J) to form a one-dimensional array, that for (3 ^> 1 has the Hamiltonian 

N N 

H = ~ Y\ J *v(°t°7-i + + Jz^U - E( A < + Bat) ■ (3) 

i=2 i=l 
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The terms J z <jfaf +l are due to the small inductive coupling between adjacent qubits. 
Here J z = 2M m I^ where M q>g is their mutual magnetic inductance. The coupling 
constant J z could in principle be increased by a common Josephson junction between 
two neighboring qubits 0. The tunneling amplitude A between the states |0) and |1) of 
the coupled qubits differs from the value A for individual non-coupled qubits, because 
coupling suppresses independent tunneling events in which only one qubit changes its 
state. Also, simultaneous tunneling events |11) « — > 1 00) for two neighboring qubits are 
suppressed and therefore we neglect such processes in our model. Correlated tunneling 
events |10) < — > |01) are unaffected by the coupling. 

The Hamiltonian (JHl) contains a term A^i^f, i-e. it does not conserve the z- 
component of the total spin (which is equivalent to the number of sites in the excited 
state |1)). Therefore, the theory proposed in pQ is not valid in our case. However, if 
(3 ^> 1 A is much less than J xy PJ and we can neglect this term at first. Later we will 
use perturbation theory to analyze how nonzero values of A affect the results. 

We assume that the gate capacitances are equal to jC. As was shown in [TT] using 
the quasiclassical approach, A can be obtained as 



Ao = y EjE c * 

where 



2(4a 2 — 1) , / — : . / r~~ arccos(a/2) , . . , 

; expi-AJMaEjUl - a 2 /4 v 1 ; ))(4) 



\ a(l+ 7 ) ^ -V— « /- 2a 



(5) 

We now want to consider two interacting qubits that are coupled by a capacitor 
PC, see Fig. |21 The dynamics of the qubit can be described by the motion of a fictitious 
particle in a potential (see Fig. HJ) with two local minima, that correspond to the states 
|0) and |1) [10J. Therefore the collective dynamics of the two qubits can be described 
by the motion of a particle in a two-dimensional potential with four minima. Without 
interaction the effective mass of the particle is M, see Eq. (|5|). When two qubits are 
connected by the capacitor, the effective mass to move in (0,0) <->• (1,1) direction is 
M + 2m* (here m* = (h/2e) 2 (3C), the effective mass for independent qubit tunneling 
events is M + m* and the effective mass for tunneling in (1, 0) <-> (0, 1) direction is equal 
to M. From these formulas one can see that the tunneling is suppressed in all directions 
except (1,0) <-> (0, 1), if m* 3> M. Due to this fact state transfer with high fidelity is 
possible. 

Using the WKB- approach and realistic qubit parameters from ^T] and ^U], namely 
a = 0.75, 7 = 0.02, we calculate A and J xy for our Hamiltonian: 



A = A exp(-0.49 V / E J /E c ( V / l + (3/5 - 1)) , (6) 



U xy = A e- a49 V^7^(l - e -o.98V^/s (Vi+^5-i)) . ( 7 ) 
With Ej/E c w 100, we obtain 



A/A = exp(-4.9( V / l + (3/5 - 1)) . (8) 
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Therefore, independent tunneling is effectively suppressed for f3 ~ 10. A and 4J xy 
coincide for (3 = 15. For (3 = 20, AJ xy is three times larger and for /3 — 30 it is 25 times 
larger than A. In this case, as we will show later, A can be neglected. 

For A = 0, the Hamiltonian Q is that of an asymmetric (XXZ) Heisenberg model 
in the presence of a magnetic field, 

N N 

Hl = ~ EI-^WVVi + vTati) + J^lx) ~ E B < ■ ( 9 ) 

1=2 8=1 

We now want to calculate the fidelity of the state transfer. The chain is initialized 
in the state |00...00) by first choosing a large negative value for the parameter B, see 
Eqs. and (J2J). Then, the first qubit is prepared in the state \ipi n ) i.e., the total state 
of the array is \tpi n , 00. ..00). This is not an eigenstate of the Hamiltonian (|9]L therefore 
the system will evolve in time. After a time t the state of the last qubit is read out. In 
general the last qubit will be in a mixed state, which is described by a density matrix 
p out . Following [T], we average the fidelity over all pure input states on the Bloch sphere 

F(t) = ±-j (ii> in \ Pmlt (tMn)dn (io) 

to obtain a quantity 1/2 < F(t) < 1 that measures the quality of transmission 
independent of \ipi n ). 



3. Calculation of the average fidelity 

We perform our calculations in the basis \k) = |00...010...0) for which the spin in the 
k-th qubit is in the state |1) and all others are in the state |0). The Hamiltonian Q 
of the array commutes with the z-component of the total spin J2i Therefore we can 
use the results of [T] to calculate the average fidelity in terms of f^ N {t) = (l\e~ lHLt \N) 
i.e., the transition amplitude of the excitation over the array. The average fidelity can 
then be expressed as 

mA 1 , l/&r(*)| a , !/&(*)! cos(7) nu 
F{t) = 2 + —^ + 3^ ' 

where 7 = Aig(ff N (t)) is the argument of the complex quantity f^ N (t). 

Varying the magnetic field one can make 7 a multiple of 2tt to maximize the 
average fidelity, such that the maximum fidelity will correspond to the maximum of 
|/ijv(t)|- Furthermore, the fidelity of any state transfer is unity, if the modulus of the 
amplitude to transmit the state |1) across the array is unity. The fidelity for a given 
state |V>m) = cos(0/2)|O) + sin(0/2)|l) in the case fg N (t) = f = \ f\ is 

F{0, <p) = ^ + cos(0)i^ + cos 2 (0)^- . (12) 

It changes monotonically from 1 for the |0) state to f 2 for the |1) state. For 
/ 7^ \f\ the fidelity can have a local minimum for 9 = arccos( ^Mjzrpn ) if 



l/l > cos( 7 )/3 + yW(7) + 3/3. 
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Figure 3. First fidelity maximum for an array with a — 0.75, 7 = 0.02, Ej/Ec = 100, 
(3 = 30 and Ej = 3GHz, a = 0.1. 
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Figure 4. Time (in units of 1/Ej) at which the first fidelity maximum is achieved. It 
is proportional to the length of the chain and depends on the coefficient J xy . 

We will now calculate l/iVWI i n ^ ne case A = 0. The eigenfunctions of Hl can be 
described as follows: 



N 



\k) = J2 b k,n\ n ) ■ 



(13) 



n=l 



From the Schrodinger equation 



H L \k) = (B(N - 2) - J Z {N - 5))|fc) - 2J Z (6 M |1) + b k>N \N)) 



N-l 



-4J xy (b kt2 \l) + J2 (Kn-l + h,n+l)\n) + 6jfc,jV-l|A0) 



(14) 



n=2 
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Figure 5. Fidelity as a function of time (in units of 1/Ej) for a chain with N = 8. 
Upper panel: first fidelity maximum at small times. Lower panel: fidelity maxima 
around t — 198. The parameters are chosen as in Fig. [3J 



we obtain the following system of equations for the coefficients bk, n 
b k , n -i + h,n+i = Dh,n (n G [2, N - 1]) 

ab k ,i + b kt2 = Db ks i (15) 
obk,N + bk,N-i = Dbk,N 

where a = J z /2J xy and D is a constant. From the first two equations b ky i can be 
expressed in terms of b^i as 

h,i = Pi(D k )b kA , (16) 
here D k , k = 1, N are the roots of 

(D-a)P N (D) =P N - 1 (D) . (17) 
Pi(D) is a polynomial, that is determined recursively 

= P 2 = D-a, Pi = - P^, i = 3,...,N. (18) 

The coefficient 5^1 can be found from the normalization conditions 

m = 4 - " L = wm + z + mD>) ■ (19) 
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Figure 6. Fidelity maxima for times less than 4000/ J z , all the chain parameters are 
as in Fig. |3 
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Figure 7. Times (in units of l/Ej) at which the fidelity maxima in Fig. are 
achieved. 



Thus we have determined the eigenfunctions of the Hamiltonian and can find its 
eigenenergies 

E k = -J z (N-5)+B(N-2)-AD k J xy . (20) 
Setting E Q = 0, we obtain 

E k = 2B + AJ z -AD k J xy . (21) 
The transition amplitude of the excitation through the array is given by 

N N 

/&(*) = Y,Cm(N\~k)e-^ = J2 b Kl b KN e-^ . (22) 

k=l k=l 

Using these formulas we have numerically calculated the average fidelities for 
different chain lengths and ratios a = J z /2J xy . The most relevant quantities for practical 
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purposes are the first fidelity maxima, see Fig. Eland Fig. that we will call "fidelity" 
in the rest of this article. 

For short-length chains the average fidelity is higher than 0.9. This makes persistent 
qubit arrays good candidates for transmission lines in quantum computers, that are 
based on flux degrees of freedom. Also they can be effectively used in the two-chain 
method proposed for achieving perfect state transfer The fidelity has a complicated 
oscillating behavior as a function of time, see Fig. El There are many local maxima, and 
the first of them is usually not the global maximum. Therefore, waiting long enough, 
we can achieve a higher fidelity. This can be seen by comparing Fig. Eland Fig. El 

However, the waiting times, i.e. the times, at which the maximum peaks of the 
fidelity shown in Fig. El occur are much longer than for the first maximum, see Fig. 
13 Therefore, from a practical point of view the first maxima in the fidelity are more 
relevant. 

Decoherence is another important reason why practical realizations of our proposal 
would have to focus on the first fidelity maximum. Like any physical realization of a 
qubit, flux qubits are characterized by a finite dephasing time, and in a recent experiment 
times of order ~ 20ns were reported for a single flux qubit [12,. Since the time for the 
appearance of the first fidelity maximum is of order hL/Ej. As a simple estimate of the 
effects of decoherence, we compare this time with the dephasing time, which leads to a a 
limit of the length of the array of L ~ r^Ej/h ~ 100. Additional maxima after the first 
one will be further reduced by decoherence since they correspond to states traversing 
the array more than once. 

To maximize the fidelity 7 = Arg(f^ N (t)) has to be chosen equal to zero. This can 
be done by varying the magnetic field, so that —2Bt + 70 = 2im. Here 70 is transition 
amplitude phase for B = 0. To achieve more control of the qubit parameters the central 
junction can be replaced by a SQUID [TT] . 

The works of Bose pQ and Christandl et al. |Hj correspond to spin chains with a 
particular form of the Hamiltonian Hi {J z /2Jxy — 1, Jz/^Jxy — 0). We have checked 
that in these limits our results agree with jT] and 

As mentioned above, the Hamiltonian of the real chain contains a term A J^Li erf, 
that does not conserve the z-component of the total spin (i.e. the number of excitations). 
A is small, however, i.e. we can use perturbation theory to analyze the influence of this 
term on the average fidelity. In this case we need to do calculations in a larger (2^ + 1)- 
dimensional space, because in principle any number of excitations is possible. One 
can easily show, that in zero-order approximation the fidelity and the N + 1 lowest 
eigenstates will be the same as in the unperturbed case. The first-order corrections are 
zero, because (k\af\k) = for the lowest eigenstates. So only the second-order terms, 
which are proportional to A 2 , affect the fidelity. The influence of the symmetry-breaking 
term therefore vanishes quadratically with A. 

From Fig. EJone can see that for qubits with the parameters mentioned in Fig. EJit is 
sufficient to choose coupling capacitors with about 25-30 times the junction capacitance. 
In this case we can neglect the influence of A. For long times this term becomes more 
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Figure 8. Fidelity dependence on (3 for chain length N = 5. 



important. This is another reason why only the first maxima are useful for practical 
realizations of high-fidelity transmission lines. One can, in principle, raise (3 to make 
the A-term less important for the maxima that occur later, but in this case the charging 
energy will increase and this will influence the fidelity and the properties of the qubit. 



4. Conclusions 



We have shown that a persistent-current qubit array is a good candidate for quantum 
state transfer with high fidelity in flux-qubit based quantum computers. For short- 
length chains the average fidelity of state transfer is higher than 0.9. Therefore, this 
type of array can be effectively used in the two-chain algorithm [Oj for achieving perfect 
state transfer. The influence of the term proportional to Aa x , that does not commute 
with the z-component of the total spin, is quadratic in A and can be neglected at small 
times for (3 ^> 1. 

We would like to thank R. Fazio for valuable discussions. 
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